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We compute fermionic observables relevant to the study of chiral sym- 
metry in quenched QCD using the Overlap-Dirac operator for a wide range 
of the fermion mass. We use analytical results to disentangle the contribu- 



tion from exact zero modes and simplify our numerical computations. Details 

ON ' 



concerning the numerical implementation of the Overlap-Dirac operator are 



p 1 presented. 
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I. INTRODUCTION 

The Overlap-Dirac operator Q derived from the overlap formalism for chiral fermions on 
the lattice @ provides a non-perturbative regularization of chiral fermions coupled vectori- 
ally to a gauge field. The massless Overlap-Dirac operator is given by (the lattice spacing 
is set to one) 

£>(0) = i[l + 75#«]- (1) 

H a is a hermitian operator that depends on the background gauge field and has eigenvalues 
equal to ±1. There are two sets of choices for H a . The simplest example in the first set is 
e(H w ) where 'j 5 H w (—m) (m c < m < 2) is the usual Wilson-Dirac operator on the lattice 



with a negative fermion mass and — m c is the critical mass, which goes to zero in the 
continuum limit. One could replace H w by any improvements of the Wilson-Dirac operator 
with the mass in the supercritical region but smaller than where the doublers become light. 
The simplest example in the other set is obtained from the domain wall formalism || and 
H a = e(\og(T w )) where T w (m) is the transfer matrix associated with the propagation of 
the four dimensional Wilson fermion in the extra direction Here m c < m < 2 is the 
domain wall mass. Again, one can replace the five dimensional Wilson-Dirac fermion by any 
improvements. The massive Overlap-Dirac operator is given by 

DM = i[l + + (l-/i) 75 # a ] (2) 

with < fi < 1 describing fermions with positive mass all the way from zero to infinity. 
For small /i, \x is proportional to the fermion mass. When H a = e(log(T M) )), the bare mass 
fi is exactly the standard mass term in the domain wall formalism HQ. The analytical 
expressions derived in this paper are valid for any generic choice of H a . Our numerical 
results are obtained with the choice of H a = e(H w ). 

The massless Overlap-Dirac operator has exact zero eigenvalues, due to topology, in 
topologically non-trivial gauge fields, and these zero eigenvalues are always paired with 
eigenvalues exactly equal to unity fl]f|]. All exact zero eigenvalues are chiral, and their 
partners at unity are also chiral but with opposite chirality. The rest of the eigenvalues of 
the Overlap-Dirac operator come in complex conjugate pairs and one eigenvector is obtained 
by the operation of 75 on the other eigenvector. The behavior of the Overlap-Dirac operator 
on the lattice is essentially that of a continuum massless fermion. 

The external fermion propagator is given by 

D-\^) = (1 - /i)- 1 [D-\v) - l] . (3) 

The subtraction at fi = is evident from the original overlap formalism |2[] and the mass- 
less propagator anti-commutes with 75 []I],|7|]. With our choice of subtraction and overall 
normalization the propagator satisfies the relation 



^\[~ f5 D-\ri]\b) = (tf\D- 1 ^)\b) Wb satisfying 75 |6) = ±|6> (4) 

for all values of fi in an arbitrary gauge field background. This relation will be evident 
from the discussion in section |T|. The fermion propagator on the lattice is related to the 
continuum propagator by 

D~ l (m q ) = Z^ l D~ x {ii) with m q = Z~ X fi (5) 

where Z m and Z^ are the mass and wavefunction renormalizations, respectively Requiring 
that (fD hold in the continuum, results in Z^Z m = 1. For H a = e(H w (m)) we find 

Z f = Z- 1 = 2m (6) 

at tree level, and a tadpole improved estimate gives 

^ = ^ l = ^-[m-4(l-«o)] , (7) 

where Uq is one's favorite choice for the tadpole link value. Most consistently, for the above 
relation, it is obtained from m c , the critical mass of usual Wilson fermion spectroscopy. 

In this paper, we study QCD in the quenched limit on the lattice using the Overlap-Dirac 
operator. Our aim is to study the behavior of chiral susceptibilities on a finite lattice in the 
massless limit. In section H, we derive analytic expressions for various chiral susceptibilities 
using the spectral decomposition of the Overlap-Dirac operator. These formulae do not 
reveal any new physics not already known from formal arguments in the continuum; however, 
they are valid for any lattice spacing, not only in the continuum limit. These formulae will 
help us simplify the numerical computations and disentangle contributions due to global 
topology from the rest of the contributions. We also discuss the massless limit of the chiral 
susceptibilities in quenched QCD. In section [TTT], we present some details of the algorithm 
used to numerically compute the chiral susceptibilities on a finite lattice. Our numerical 
results are presented in section [TV] and compared with the analytical expressions obtained 
in section |IJ. Similar numerical work has been done recently using domain wall fermions 
with a finite extent in the fifth direction 101. 



II. CHIRAL CONDENSATE AND SUSCEPTIBILITIES 



An analysis of the Overlap-Dirac operator can be performed in an arbitrary gauge field 
background by noting that H a can be written in block diagonal form with 2x2 blocks in 
the chiral basis @. We start by recalling a few properties of the Overlap-Dirac operator 
D for a single Dirac fermion: ( / y 5 D) 2 commutes with 75 and both operators can therefore 
be diagonalized simultaneously The massive Overlap-Dirac operator can be related to the 
massless one, from (E3), as 



D(p) = (l-p) 



D(Q)+ 



(8) 



We find it convenient to work in the chiral eigenbasis of (7s-D(0)) 2 for massless fermions, 
jj, = 0. In this basis D(fi) takes on the following block diagonal form ||: 

• There are \Q\ blocks of the form 

> 0' 

1 

associated with the global topology of the configuration. These \Q\ blocks will be 
robust under perturbations of the gauge field background. The two modes per block 
have opposite chiralities. The chiralities of the modes with eigenvalues \i will all be 
the same and dictated by the sign of the global topology Q. 

There are 2NV — \Q\ blocks of the form 



;i-^)a?+// (i-//)Wi-A 



-(1 - n)\\Jl — Xj (l-At)A? + /i 
with < Aj < 1, i = 1, 2, • • • , 2NV — \Q\, depending on the background gauge field. 
The Af's are the eigenvalues of (7 5 D(0)) 2 . We leave open the possibility that n' z of 
these eigenvalues could be exactly zero. These are accidental zeros in that they would 
be lifted by a small perturbation of the gauge field background. In each subspace, 

75 ^(o -1 
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The matrix can be inverted in the above form to give the following block diagonal form for 
the external propagator = (1 — /i) -1 [Z) _1 (/x) — 1]: 

• There are \Q\ blocks of the form 

i 0' 


• There are 2NV — \Q\ blocks of the form 

i Mi -a?) -kji^xi 

We will use (• • -)a to denote the expectation value of a fermionic observable in a fixed 
gauge field background, and (• • •) for the expectation value averaged over a gauge field 
ensemble. From the above expression for the fermion propagator in a background gauge 
field, we obtain 

\ E<^)75^)>a = ^Tr^D- 1 ] = -5- ( 9 ) 

in a gauge field background with global topology Q. Here Tr denotes a trace over space, color 
and Dirac indices, (^j^ip) goes to zero when averaged over all gauge field configurations. 
We also find 

i WW* - = $ + 7 ? MT^h ■ (10) 

To discuss the gauge field average of (ipip) in the massless limit, we have to distinguish four 
cases []: 

• We consider full QCD with n/ > 1 flavors in a fixed finite volume. Gauge fields with 
non-zero topology are then suppressed by jx n ^\ the contribution from the zeromodes 
to det D(fi). Accidental zeromodes are also suppressed by the fermion determinant, 
by fi 2n f n ' z . Therefore, we find from eq. (|10|) 



1 We do not consider the case of one flavor QCD since the chiral symmetry is broken by the 
anomaly. 
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±J2$(xW(x)) = (11) 

' X 

and the chiral condensate vanishes in the massless limit in a finite volume. 

We consider full QCD with n/ > 1 flavors in the infinite volume limit, taking the 
volume to infinity before taking the mass to zero. In this limit we expect (\Q\) to 
be proportional to \/V, and therefore the first term in fln]) will go to zero in the 
infinite volume limit. The spectrum, X i} becomes continuous in the infinite volume 
limit, represented by a spectral density, p(A). Then we obtain, 

lim lim I]T^(a;)V>(2)> = tt P (0) (12) 

and the chiral condensate is proportional to the spectral density at zero just as expected 
in the continuum ||. 

We consider a quenched theory in a fixed finite volume. Now, gauge field configurations 
with nontrivial topology and accidental zeromodes are not suppressed by a fermion 
determinant, and we obtain 

lim - yOip(xWx)) = lim + 2 <) (13) 

2n' z is the number of accidental, paired zero modes. Since they are lifted by a small 
perturbation of the gauge fields, the gauge fields for which such accidental zero modes 
exist are probably of measure zero, and we expect (n' z ) = 0. Thus in any finite volume 
(ipip) diverges like - in the quenched theory with a coefficient equal to ^M-. The 
expectation value here is over a pure gauge ensemble, and the coefficient is therefore 
expected to be proportional to l/y/V. 

We consider a quenched theory in the infinite volume limit, taking the volume to 
infinity before taking the mass to zero. Then, the first term in (|T0| ) vanishes and from 
the second we obtain 
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lim lim i Y,^( x )^( x )) = *P( Q ) ( 14 ) 

where p(0) is the spectral density at zero in a pure gauge ensemble. 

A two flavor theory has a flavor SU(2)yX SU(2)^ =0(4) symmetry. Consider the (real) 
0(4) vector made up of 



= 0" = ipiix^r^jix), f = Y, M x )M x )) ( 15 ) 

ij i 

where fo is also known as a, and the opposite parity (real) 0(4) vector made up of 

4> = K = ~Yl^ x ) T ij^A x )i v = iE^W^W) • ( 16 ) 

ij i 

Here r a , a = 1,2,3, are the SU(2) generators. All eight components are invariant under a 
global U(l) rotation associated with fermion number. Under a vector SU(2) flavor rotation, 
the ir a rotate among themselves, the Oq rotate among themselves, and /q and t] are left 
invariant. Under an axial SU(2) flavor rotation, the 7t a mix with /o and the mix with r\. 
Finally, under a global U(l) axial rotation, mixes with 0. 
Consider the space averaged pion propagator 

X, = E ^K(*K(2/)> = |(Tr(75^)- 2 ) • (17) 

x,y V ' 

(7 5 _D)~ 2 is diagonal in the chiral basis. There are \Q\ eigenvalues equal to -\ and \Q\ zero 
eigenvalues from the global topology. The non-zero eigenvalues appear in multiples of two 
and are equal to (1 — A 2 )/ (A?(l — jj 2 ) + jj 2 ) which is nothing but the diagonal entries of 
-D~ l . Relation (f|) in section | therefore follows. We also obtain 

Ity^ d)- 2 = M + I v 2(1 " A?) 
V 175 ' " /i 2 V v y \ 2 (i - 2 ) + /" 2 

= -(W>)a • (18) 

Therefore, we find the relation 

x, = Y.U* a ( x > a (y)) = -^)- (is) 



We note that this relation holds for any volume and /i, configuration by configuration. In 
fact, from (|J) it follows that the relation holds for any chiral random source that might be 
used for stochastic estimates of Xn and (ipip). 
Using, similarly, 

\r(n^ IGI 1 v 2(1 - Ag - A?) - A?] 

^ ) = ^ + U? [A?(l-^)+^] 2 (20) 

we find for the space averaged ao propagator 

X«o = E ^<oS(*)aS(l/)> = -|(M- 2 ) = 2{Uw)a) ■ (21) 
For the space averaged /o propagator we find 

Xf = E ^(Mx)Mv)) = ^D- 1 ] 2 ) - ^D- 2 } , (22) 

X,j/ * V V 

and for the space averaged 77 propagator 

X, = E^(v(xHy)) = f (Tr^r 2 ) - ^([Tr^)- 1 ] 2 ) . (23) 

To discuss these chiral susceptibilities, we again distinguish four cases: 

• We consider full QCD in a finite volume and recall that gauge fields with non-zero 
topology are suppressed by yU 2 ' Q ', while gauge fields with accidental zeromodes are 
suppressed by /i 4n ' z from the fermion determinant. Using eqs. fll8D , (p0|) and (H) we 
find 

. 8/ A» 2 (l - Ag) 2 \ 4(|Q1) 4(Q 2 ) 

X, X„, v\y[a?(i-/< 2 )+/< 2 ] 2 / M 2 f M 2 V ' 

Taking the /i — > limit, the first term on the right hand side vanishes, while for the 
others only the sectors with Q — ±1 contribute, for which the two remaining terms 
cancel. Therefore we obtain 

lim , iXv ~ Xao) = , (25) 
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and, as expected in a finite volume, the 0(4) symmetry is unbroken in the massless 
limit. On the other hand, consider 

2 d - 

U = Xn~ Xa = -(lH>) - 2(-T-(W)a) 

\i dfj, 
8/ A« a (l-\ ? ) a \ , 4(131) 

= X,-Xa + ^/- ■ (26) 



We have used (|T9|) and (pi) in the first line, (|i8|) and (|20|) in the second line and (I 



in the last line. Since Xv ~ Xa vanishes in the massless limit, we obtain 

4(Q 2 ) 



limw = lim^/. (27) 



The U(1)a symmetry, therefore, remains broken in a finite volume, due to topology. 
Only the sectors with Q = ±1 contribute in the massless limit, giving a finite result. 
In fact, the four fermi operator making up Xn ~ Xa written explicitly in terms of left 
and right handed spinors, is nothing but the t'Hooft vertex. 

We consider full QCD in the infinite volume limit, taking the volume to infinity first. 
Then, one expects the 0(4) symmetry to be broken down to 0(3)=SU(2), where the 
pions are the Goldstone bosons associated with the spontaneous symmetry breaking. 
Eq. (|T9"D is consistent with this expectation, implying that oc fi. The rj, on the 
other hand, is expected to remain massive due to the axial U(l) anomaly. Therefore 
we find, using fl23|), (|18]) and (|), 

= lim lim fi Xv = hm lim ^^(^Mv)) 
= lim lim ( 2 <#>-^U , (28) 

which leads to the relation 

2(Q 2 ) 

lim lim = lim lim (ipip) > (29) 



a relation that has been derived before in the continuum 10 



We consider a quenched theory in a fixed finite volume, assuming two flavors of valence 
fermions. Now, gauge fields with non-zero topology are not suppressed, and we obtain 



from (f24j) for x-q ~ Xa in the massless limit: 



lim ( Xv - Xa ) = lim -^{\Q\ + 2< - Q 2 ) . (30) 

Though we expect (n' z ) = for the number of accidental paired zero modes, we find 
not only that the 0(4) symmetry remains broken in the massless quenched theory in 
a finite volume, but that Xv ~ Xa diverges in the massless limit, unless we restrict the 
quenched theory to the gauge field sector that has \Q\ = 0, 1. With this restriction, 
one can still see the global U(l) anomaly; however, u diverges as 4(|Q|)/(/i 2 V r ) in the 
massless limit. 

We consider a quenched theory in the infinite volume limit, taking the volume to 
infinity before taking the mass to zero. From (|I7]) and ([19]) it is clear that 



Xn-Xv = (^7> ■ (31) 



If we restrict ourselves to the topologically trivial sector 7r and 77 will be degenerate. If 
we allow all topological sectors, the right hand side of the above equation will diverge in 
the massless limit and the rj particle will not be well defined (x v will become negative). 

III. NUMERICAL DETAILS 

The main quantity that needs to be computed numerically is the fermion propagator 
D~ l . Certain properties of the Overlap-Dirac operator enable us to compute the propagator 



for several fermion masses at one time using the multiple Krylov space solver |TT] and also 
go directly to the massless limit. Our numerical procedure to compute the fermion bilinear 
in (pTJ| ) and u in ( p6|) proceeds as follows. 

• We note that 
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,2 



1 — n z 



(32) 



with 

1 1 



H 2 o (0) = ^ + l [l 5 e(H w ) + e(H w ) l5 ] (33) 



Eq. ( P2] ) implies that we can solve the set of equations H 2 (fi)rj(fi) = b for several 
masses, /x, simultaneously (for the same right hand b) using the multiple Krylov 



space solver described in Ref. ||11|| . We will refer to this as the outer conjugate 
gradient inversion. 

— It is now easy to see that [H^(ji), 75] = 0, implying that one can work with the 
source b and solutions 77 (/i) restricted to one chiral sector. 

The numerically expensive part of the Overlap-Dirac operator is the action of H^(0) 
on a vector since it involves the action of \y^,e{H w ) + e{H w )^f^\ on a vector. If the vector 
b is chiral {i.e. 756 = ±6) then, [y 5 e(H w ) + e(if„,)7 5 ]6 = [75 ± l]e(H w )b. Therefore we 
only need to compute the action of e(H w ) on a single vector. 

In order to compute the action of e(H w ) on a vector one needs a representation of the 
operator e(H w ). For this purpose we use the optimal rational approximation for e(H w ) 
described in Ref. ||. As described in 0, the action of e(H w ) in the optimal rational 
approximation amounts to solving equations of the form (H^ + qk)rjk = b for several 



k which again can be done efficiently using the multiple Krylov solver |TTJ. We will 
refer to this as the inner conjugate gradient inversion. 

A further improvement can be made by projecting out a few low lying eigenvectors of 
H w before the action of e(H w ). The low lying eigenvectors can be computed using the 
Ritz functional |T^]. This effectively reduces the condition number of H w and speeds 
up the solution of (H 2 +q k )rj k = b. In addition the low lying subspace of H w is treated 
exactly in the action of e(H w ) thereby improving the approximation used for e(H w ). 



• Since in practice it is the action of H^(fi) on a vector we need, we can check for the 
convergence of the complete operator at each inner iteration of e(H w ). This saves some 
small amount of work at \i — and more and more as [i increases, while at /x = 1 
(corresponding to infinitely heavy fermions) no work at all is required. 

We stochastically estimate TrD -1 , Tr(7s.D) _2 and TrD~ 2 in a fixed gauge field back- 
ground. To be able to do this efficiently, and in order to compute this at arbitrarily small 
fermion masses, we first compute the low lying spectrum of 75_D(0) using the Ritz func- 
tional method [|l^] . This gives us, in particular, information about the number of zero 
modes and their chiralityf] In gauge fields with zero modes we always find that all \Q\ 
zero modes have the same chirality. We have not found any accidental zero mode pairs 
with opposite chiralities. We then stochastically estimate the three traces above in the 
chiral sector that does not have any zero modes. In this sector, the propagator with any 
source is non-singular even for zero fermion mass. Given a Gaussian random source b with 
a definite chirality we solve the equation H (n) 2 r](fi) = b for several values of \x. Then we 
compute (1 — = (fi)rj(fi) — b. The stochastic estimate of TyD^ 1 is obtained by 

computing b'£([j,) and averaging over several independent Gaussian chiral vectors b. The 
stochastic estimate of Tr(7 5 D)~ 2 is similarly obtained by computing £t (//)£(//) and averag- 
ing over several independent Gaussian chiral vectors b. To estimate Tr.D~ 2 we also compute 
(1 — fi)^'(fi) = D(fi)rj(fi) —b. Then TrD~ 2 is stochastically obtained by computing £'*(/x)£(/i) 
and averaging over several independent Gaussian chiral vectors b. Since the stochastic esti- 
mates were done in the chiral sector that does not contain any zero modes, the result has to 
be doubled to account for the contributions from the non-zero modes in both chiral sectors. 
The contributions from the zero modes due to topology, finally, are added on analytically 
from the appropriate equations in section [TT]. Since modes with A 2 = 1 do not contribute at 



2 The topology and chirality information can also be obtained from the spectral flow method of 
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all to D \ we do not have to worry about double counting these edge modes. 

IV. NUMERICAL RESULTS 

For our numerical simulations, we used the usual hermitian Wilson-Dirac operator 
H(m) = ^^W^—m) with a negative fermion mass term JIJ at m — 1.65. Our main re- 
sults are for quenched SU(3) ft = 5.85, where we have three volumes, and ft = 5.7, where 
we considered one. In Table |, we show the simulation parameters and number of configu- 
rations used. From our previous work with the spectral flow of the hermitian Wilson-Dirac 
operator flnf| , we have higher statistics estimates of (\Q\) and (Q 2 ) which we quote in Table | 



with N top configurations used. For our stochastic estimates of the various traces, we used 
-Wconf configurations and 5 stochastic Gaussian random sources. We found 5 sources to be 
sufficient. For the plots of (ipip) / \i and u, we used iV con f configurations for the topology term 
with all statistical errors computed using a jackknife procedure. 

The fermion masses we used for our simulations range from to 0.5. We note that the 
number of (outer) conjugate gradient iterations we require is typically about 50, 90, and 
200 with a convergence accuracy of 10~ 6 (normalized by the source) for our three volumes 
(in increasing order) at ft = 5.85. This is consistent with our expectation that while we 
have chosen a chiral source that has no overlap with the zero modes, the smallest non-zero 
eigenvalues of H 2 (0) are decreasing with increasing volume. The number of inner conjugate 
gradient iterations (used for the application of e(H w ) on a vector) is about 100, 105, and 150 
with a convergence accuracy of 10~ 6 for our three volumes at ft — 5.85. We projected out the 
10, 10, and 20 lowest eigenvectors of H w in the computation of e(H w ), respectively, before 
using the optimal rational approximation to compute e(H w ) in the orthogonal subspace. 



While costly, the pole method ||[T5| of implementing e(H w ) is clearly the most cost effective 
in terms of floating point operations compared to other implementations of the Overlap- 
Dirac operator P.|14||. 



We show our result for {ipip) without the topology term included in Figure [l]. Shown 
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(but not discernible) are the results for all four volumes. An expanded view for small mass 
is shown in Figure 0. For small fi, a linear dependence on \i with zero intercept is expected 
from Eq. (|i~0|). If there is to be a chiral condensate in the infinite volume limit the slope 
near fi = should increase with volume. To see this effect we show (ipip) / /i for all our data 
sets in Figure |3|. For sufficiently small /i, we see linearity in /i due to finite volume effects. 
This linearity is manifested as a constancy in {ijjip) / [i. As the volume is increased, the mass 
region where linearity sets in is shifted to smaller values. However, the slope near /i = does 
not increase proportional to the volume and we conclude that there is no definite evidence 
for a chiral condensate in the infinite volume limit yet. 

We show (ipip) with the topology term included in Figure f|. As expected from Eq. ( |10| ) . 
the topology term diverges in a finite volume for fi — * 0. However, the singularity is decreas- 
ing for increasing volume. For sufficiently large volumes, we expect that the quenched 
scales like \/V . From Table | we see that the increase of (\Q\) from 8 3 x 16 to 8 x 16 3 is 
consistent with a \/V growth. But this is not the case when we compare the value at 6 3 x 12 
with the one at 8 3 x 16. This is attributed to a finite volume effect present at 6 3 x 12 

More insight into the possible onset of spontaneous chiral symmetry breaking can be 
obtained by studying uo in Eq. (^). If chiral symmetry is broken then u should diverge 
as -; otherwise, it should go to zero as fi 2 . Evidence for chiral symmetry breaking in the 
infinite volume limit should show up as a - behavior in u at some moderately small fermion 
masses ( an increase in the value of u as one decreases the fermion mass) that turns over to 
a fi 2 behavior at very small fermion masses. We show u without the topology term added 
in Figure |5|. We see a turnover developing in the j3 = 5.85 data as the volume is increased 
around fi ~ 10~ 2 . The turnover is also present in the data at (3 — 5.7. Therefore it is possible 
that there is an onset of chiral symmetry breaking, but the effect is not strong enough for 
us to be able to extract a value for the chiral condensate. 

When we add the topology term in a;, as shown in Figure |6|, the turnover region is com- 
pletely obscured. This result demonstrates how finite volume effects relevant for a study of 
chiral symmetry are obscured by topology. To be able to extract the chiral condensate in a 
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quenched theory from finite volume studies it would therefore be helpful to remove the con- 
tribution from topology. This is possible with the Overlap-Dirac operator as demonstrated 
here. 

V. CONCLUSIONS 

We have derived the standard continuum relations among various fermionic observables 
by working on a finite volume lattice with the Overlap-Dirac operator. Our results are 
general and apply equally well to the Overlap-Dirac operator of Neuberger [0 and domain 
wall fermions for infinite extent in the extra fifth dimension Q. These relations do not reveal 
any new physics but allow us to disentangle the contribution of topology and simplify our 
numerical simulations. We should emphasize that the relations follow only because Overlap 
fermions satisfy the usual chiral symmetries at any lattice spacing, not only in the continuum 
limit. 

We have shown that the standard relation between the pion susceptibility and the fermion 
bilinear (c.f. ([H])) is properly reproduced by our implementation of the Overlap-Dirac oper- 
ator even in the chiral limit. We found that for the volumes studied the chiral condensate in 
the quenched approximation at small fermion mass is dominated by the contribution from 
the global topology. While this contribution eventually vanishes in the infinite volume limit, 
it overpowers any signs for an onset of a non- vanishing condensate for the volumes that we 
were able to study. Removing this topological contribution, which is possible in our imple- 
mentation of the Overlap-Dirac operator, enabled us to perform a finite volume analysis of 
chiral symmetry breaking. Still, our numerical simulations for quenched SU(3) gauge theory 
do not present strong evidence for a chiral condensate in the infinite volume limit of the 
quenched theory; however, there is some evidence for the onset of chiral symmetry breaking 
as the volume is increased. 
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size 




iVtop 


(IQI> 


(Q 2 ) 


5.85 


8 x 16 3 


31 


40 


2.54 ±0.35 


10.7 ±3.5 


5.85 


8 3 x 16 


50 


200 


1.115 ±0.067 


2.13 ±0.22 


5.85 


6 3 x 12 


97 


400 


0.268 ± 0.025 


0.323 ± 0.037 


5.70 


8 3 x 16 


50 


50 


1.98 ±0.23 


6.6 ± 1.3 



TABLE I. Simulation parameters at m = 1.65. The number of configurations iV con f was used 
for the trace estimates, and iVt op configurations were used for the estimates of (\Q\) and (Q 2 )- The 
spectral flow method of ref. []ll|] was used for the estimates of (\Q\) and (Q 2 ) for all but the first 
line. 
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FIG. 1. Plot of (if)ij)) without topology for the data sets in Table 



0.6 



19 



0.05 



0.00 



i 



si 



1 



CD 



O = 



55, 8x16^ 



x £ = 5.85, 8 xl6 
o (6 = 5.85, 6 3 xl2 
x 13 = 5.70, 8 3 xl6 







0.00 



0.01 



0.02 



FIG. 2. Expanded view of (V'V') without topology. 
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FIG. 3. Plot of without topology. Strong finite volume effects are manifested as con- 

stancy in jjl. The mass region where finite volume effects sets in moves to smaller mass values as 
the volume is increased. 
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FIG. 4. Plot of (V>V>> with topology. 
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FIG. 5. Plot of w without topology. 
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FIG. 6. Plot of a; with topology. 
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